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ABSTRACT 

A hydrostatic stress-dependent, anisotropic model of viscoplasticity is formulated as an extension 
of Bodner's model. This represents a further extension of the isotropic Bodner model over that made to 
anisotropy by Robinson and Miti-Kavuma (1993). Here, account is made of the inelastic deformation that 
can occur in metallic composites under hydrostatic stress. A procedure for determining the material 
parameters is identified that is virtually identical to the established characterization procedure for the 
original Bodner model. Characterization can be achieved using longitudinal / transverse tensile and shear 
tests and hydrostatic stress tests; alternatively, four off-axis tensile tests can be used. Conditions for a yield 
stress minimum under off-axis tension are discussed. The model is applied to a W/Cu composite, 
characterization is made using off-axis tensile data generated at NASA Lewis Research Center ( LeRC). 


1. INTRODUCTION 

The model of viscoplasticity introduced by Bodner and Partom (1975) and 
reviewed by Bodner (1987) is one of the class of unified constitutive models. It applies to 
isotropic viscoplastic solids and is particularly effective in representing rate-sensitive, non 
isothermal responses typical of the histograms of rocket engines, the space shuttle main 
engine (SSME) being a prominent example. A major strength of Bodner's model is its 
simplicity and the relative ease with which the material parameters can be determined 
through experiment. 

The Bodner model has been extended by Robinson and Miti-Kavuma (1993) to 
represent anisotropic materials with application to metal matrix composites (MMC s). 
Their extended model is capable of representing strong initial anisotropy yet is based on a 
scalar state variable under the assertion that induced anisotropy, through directional or 
kinematic hardening, is negligible compared to the strong initial anisotropy. Robinson and 
Miti-Kavuma (1993) account for varying fiber volume fraction through the introduction of 
a nonlinear rule of mixtures. They apply the anistropic model to a SiC/Ti composite. 

It has been pointed out by Dvorak and Rao (1976), Dvorak (1983) and Rogers 
(1990), among others, that plastic deformation (and damage) in the matrix of strongly 
reinforced composites can result from the application of hydrostatic stress. This suggests 
that the stress dependence in phenomenological constitutive relations for MMC's should 
be based not solely on deviatoric stress but also on hydrostatic stress, viz, /, =tra = <y, r 
The objective of this paper is to provide a further extension of the isotropic Bodner 



theory of viscoplasticity to include not only anisotropy, as in Robinson and Miti-Kavuma 
(1993), but an appropriate dependence on hydrostatic stress. 

Several constitutive models accounting for the influence of hydrostatic stress on 
inelastic deformation have been proposed, these include Alhrot, Boehler and Sawczuk 
(1977) for anisotropic rock, Doraivelu et al. (1984) for porous powder-metallurgical 
materials, Rogers (1990) for 'constraint hardening' fiber-reinforced composites and Chan 
et al. (1992) for deformation and damage in rock salt. 

We formulate the anisotropic, hydrostatic stress-dependent model in terms of an 

effective stress that depends on /, and invariants reflecting transverse isotropy, cf., Lance 

and Robinson (1971) and Robinson and Duffy (1990), thus replacing.^L/j’ in the 
isotropic Bodner model. The multiaxial model is specialized for longitudinal / transverse 
tension and shear and hydrostatic pressure; these constitute the fundamental or natural 
stress states (or tests) for the model. Relationships between material constants specific to 
these natural states and those of the multiaxial model are identified such that the governing 
equations in each case reduce identically to the form of the uniaxial equations of the 
isotropic Bodner model. This enables the same well established characterization procedure 
developed for the original Bodner model to be applied to the present anisotropic model. A 
complete procedure for determination of the material parameters is specified on this basis. 
Thus, the extended model retains the simplicity of the Bodner theory in the ease with 
which the material constants are determined. 

An application of the model is made to a unidirectional W/Cu (Tungsten/Copper) 
composite. Characterization is based on data generated by M.J. Verilli of the NASA /Lewis 
Research Center, the testing temperature is 260C and the fiber volume fraction is » 9% . 


2. THE MULTIAXIAL MODEL 


We state the multiaxial form of the anisotropic, hydrostatic stress-dependent model 
for isothermal conditions. Temperature dependence can be included as in Bodner (1987) 
and Robinson and Miti-Kavuma (1993). The flow law is 

4 ( 1 ) ] r ‘ 

(2.1) 

and the evolution law is 

Z = tn{Z s — 

(2.2) 


in which s v are the components of inelastic strain rate, cr is the effective stress 

5 = M7, -{(J-Jl)-(C- nVl (2.3) 

and Z is the internal state variable (stress history parameter). 
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The invariants J 2 ,J 0 ,J and /, in (2.3) are defined as 


, _ 1 

•A “ 2 SySji 


Jo ~ tyjSjj 


J = DgSpSto /, -tr tt 


(2.4) 


Also, Tjj—Sjj £( An s jk A* s h 2J 0 D iJ ) 2(£ 77)^0 (Ay 

+f«'-'W s (2.5) 

with 0 £ £ < 1 0 < 77 < ^- 0<f<l (2-6) 

o-„ are the components of stress, s v are the components of deviatoric stress and 

D . = dd is a symmetric orientation tensor, cf., Spencer (1972) and Robinson and Duffy 

(1990), formed by the product of the components of a unit vector d t that designates the 
local direction of transverse isotropy. 

The quantities 

n,m,'s g ,Z, and Z 0 ( the initial value of Z) (2.7) 

are the material constants of viscoplasticity; £, 77 and Q are constants specifying the degree 
of anisotropy and the hydrostatic stress dependence. 

3. SPECIAL CASES 


Hydrostatic Stress Independence (£—4 77 — > 0) 

The model (2. l)-(2.6) reduces to one insensitive to hydrostatic stress as 
<■-477 -> 0. Then, (2.3) and (2.5) become 


and 


a= 3 


j 2 -rj 
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(3.1) 

(3.2) 


which are identical to those proposed earlier in the anisotropic extension of Bodneris 
model by Robinson and Miti-Kavuma (1993). Hydrostatic stress independence implies 
inelastic incompressibility, i.e., 
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r„ = o 


and 


e u = 0 


(3.3) 


as is readily shown from (2. 1) and (3.2). Further, > 1 in (3.1) and (3.2) corresponds to 
longitudinal inextensibility, i.e.. 


A^ = ° 


(3.4) 


Equal Behavior in Longitudinal and Transverse Shear (4 = 0) 


It is often observed that the inelastic responses in longitudinal and transverse shear 
are not significantly different in some strongly reinforced metallic composites. The limiting 

case of these responses being equal corresponds to £ = 0, giving 


a= 3 




(3.5) 


r, =s„ -\s t )+~((-4nV,s f 


(3.6) 


If, in addition, the material is insensitive to hydrostatic stress (C - 4 r; -> 0 ), we have 

(3.7) 






(3.8) 


as used earlier by Robinson and Miti-Kavuma (1993) to characterize a SiC/Ti composite 
material. Further, taking 1, i.e., longitudinal inextensibility, a model analogous to that 
discussed in Binienda and Robinson (1991) results. 

Isotropy (£ = *7 = £ = 0) 

In the limiting case of full isotropy, the present model (2.1)-(2.6) reduces to that of 
Bodner and Partom (1975) and Bodner (1987) with 

cr = V^ (3-9) 


T v =s v 


(3.10) 
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4. FUNDAMENTAL STRESS STATES: DEFINITION OF MATERIAL 

PARAMETERS 


The uniaxial form of the isotropic Bodner model is 



Z* =m(Z* -Z*)ae (4.2) 


cf. Bodner (1987). a is the uniaxial stress and e is the inelastic strain rate in the stress 
direction. 


The material constants 

n,m,e* 0 ,Z* s and Z* (4.3) 

are found by correlating calculated responses based on (4.1) and (4.2) with experimental 
data obtained at strain-rates (and temperatures) of interest. The characterization procedure 
for the Bodner model is well established and discussed in Bodner and Partom (1975) and 
Bodner (1987). 

We now apply the present model (2.1)-(2.6) to the four fundamental stress states 
shown in Fig. la- Id; (a) longitudinal tension (LN), (b) transverse tension (TN), 

(c) longitudinal shear (LS) and (d) transverse shear (TS). In each case, the flow and 
evolution laws reduce identically to the form (4. l)-(4.2) provided the following 

identifications are made relating the state variable Z and the constants £ 0 ,Z S and Z 0 in 
(4 1)_(4,2) to Z and the corresponding constants s 0 ,Z s and Z 0 in the multiaxial model 
(2.1)-(2.6): 


(a) Longitudinal tension (LN1 


Z* =Z W = 




e 0 =e?=e 0 ^ 
fhl Transverse tension (TN) 


' s.O 


V* _ 7LN _ 

~VP? 


z* = z™ = 


VPt7 


e 0 = e™ = foVl-P 


7* __ 7 77^ _ 

Zi0= l0 "VTF 


(4.4) 

(4.5) 


(4.6) 

(4.7) 
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(c) Longitudinal shear (LS) 


Z* = = 


>/3(W) 


(d) Transverse shear (TS) 


^I.o 


7* = 7 U =. 

'•° "° (Ft) 


z* = z r5 = A= 
S 


e = e TS = 
6 0 c 0 


V3. 


7 * _ 7 is _ ^5,0 

■ “ “ V3 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


Because the governing equations in each case (a)-(d) reduce identically to the Bodner 
form (4.1)-(4.2), it follows that the characterization procedure used for the isotropic 
Bodner model along with uniaxial test data can be applied to the present anisotropic model 
with test data for the stress states (a)-(d). As indicated earlier, this characterization 
procedure is well established. 

Once the appropriate experiments are conducted under the stress states (a)-(d) and 
the material constants (4.4)-(4. 1 1) determined from the known Bodner characterization 

procedure, tj and £ are then determined from 


£=1- 


( yTS V 

± 0 

7 is 

V^o J 


77=1-3 


3 

2, 


'TS \ 2 

•7V 
0 ) 


<r=i-3 


(7 TS v 

Z LN 

0 J 


(4.12) 


Another fundamental stress state is that of hydrostatic pressure (Fig. le). This , 
too, leads to governing equations identical in form to (4.1) and (4.2) with a replaced by 
the hydrostatic stress p and e replaced by e = -e„, the rate of volume change. In this case 
we have 


Z* = Z H = 


V3(<r-4)7) 


(4.13) 


bl=sZ=eJ 


7 * _ yH __ 0 


(4.14) 


leading to 


C-4tj 


'gj 

2 o H ) 


(4.15) 
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The hydrostatic test. Fig. le, can replace either of the uniaxial tests Fig. la or lb as one of 
the fundamental tests, with (4. 1 5) then replacing either of the last two equations (4.12). 

Although, in principle, the tests of Fig. la-le are the natural ones for the model 
(2.1)-(2.6), the shear tests (Fig. lc and Id) and the hydrostatic test (Fig. le) are difficult 
to conduct and such data are seldom available. In the next section we consider an 
alternative set of uniaxial tests for determining the material parameters. 

5. OFF-AXIS UNIAXIAL TESTS : CHARACTERIZATION 


We now consider uniaxial tensile stress applied at an arbitrary angle 0 with the 
axis of transverse isotropy. Ideally, such data can be generated by axially loading a thin 
walled tube as shown in Fig. 2. 

The governing equations for the response s along the tensile stress cr direction 
(axial in Fig. 2) are found from (2. l)-(2.6) and, again, they reduce to the form (4. l)-(4.2), 
viz.. 


with 


s 

1 (z e> \ 

2n~ 

~ - ex p 




4 

2 V a ) 



Z e =m(Zf-Z e )as (5.2) 

z e =-r =£ — = (5.3) 




z 9 z -‘ 

" jme.’i.o 


(5.4) 


and /(«,£ = l-3(co!? 0sm J 0+i(l7-C)(2«»’0-sil' ^ 0> , 

+|(f-4u) (55) 

As argued earlier, because the governing equations reduce to the same Bodner form, we 
can determine the material parameters (5.4) using the known characterization procedure. 
We note that 0 = 0° and 0= 90° in (5.1)-(5.5) correspond, respectively, to cases (a) 
longitudinal tension (LN) and (b) transverse tension (TN) of the previous section. 

Now we suppose data has been generated as in Fig. 2 for 0=0 °,#?,#° and 90°. 

We further suppose that Z™ ,Z™ ,Z°' and Z® 1 have been determined using the Bodner 
characterization procedure. Then, from (4.5),(4.7),(5.4) and (5.5) there results, 
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fzr ) 

2 

1-7 

(5.6) 

UrJ 


7 ts \ 
^0 

7,o 

(5.7) 

7*i 

V^o ) 

i-C 

(?LN^ 

f&.ln, 0 

(5.8) 

Uo’J 

l-<r 


providing three linear equations for the determination of the anisotropy parameters tj 
and 

The remaining viscoplasticity parameters can be obtained from either of the four 
testing orientations, e.g., 0=0 0 longitudinal tension (LN), in which case we have n,m,ef 

and Zf in addition to . Using (4.5), the parameters (2.7) are known and the 

multiaxial model (2. l)-(2.6) is fully specified. 

The characterization procedure discussed above is applied to a W/Cu metallic 
composite in the next section. First, however, we examine some features of the function 

/(*,£» 7,0 in (5.3)-(5.5). 

It is easily shown that f ( 6 , £, tj,Q achieves a maximum at 


0 m = tan-‘ 


2(<T-77)-3£ 


(5.9) 


while, for example, Z* in (5.4) achieves a minimum. As Z® can be considered roughly the 
yield stress in the stress direction 0, it follows that the yield stress is a minimum at 0 m . 

In the limiting cases 


4(C-*7)-3£ = 0 (5.10) 

2(C-*7)-3£=0 (5.11) 


0 m = 0° or 90° , respectively. For the special case £ = 0, cf., (3.5) and (3.6), 0 m « 54.73° 
for all £ and tj . With 

4(C-7)~3£>0 (5.12) 

and 2(C-77)-3£<0 (513) 

0 m is imaginary (see Fig. 3) and the yield stress decreases monotonically as 0 = 0° — > 90°, 
with no minimum in 0° < 0 < 90°. 


8 



6. APPLICATION TO W/Cu 


Application of the model (2.1)-(2.6) is made to a W/Cu (Tungsten /Copper) 
composite based on data generated by M.J.Verilli ofNASA/LeRC. The testing 
temperature is 260C and the fiber volume fraction is ~ 9%. Most of the tensile data used 
were generated at a strain rate of ~ .001/s, with additional data at ~ . 000001/s used for 
establishing the strain rate dependence. Extension to other temperatures and fiber volume 
fractions follows as in Robinson and Miti-Kavuma (1993). 

In addition to longitudinal 0° (LN) and transverse 90° (TN) tensile tests, off-axis 
tensile tests with 0=5° and 15° were used as part of the data base. The off-axis tests were 
made not on thin-walled tubes as illustrated in Fig. 2 but instead on laminated cross-ply 
coupon specimens. We assume that inter laminar shear at these relatively small cross-ply 
angles has only a second order effect on the off-axis tensile response. 

Application of the characterization procedure described in the previous section 
results in the material parameters characterizing the 9% W/Cu composite at 260C shown 
in Table 1. 

Correlation of tensile data and calculations based on the material parameters of 
Table 1 is shown in Fig.4. Data and calculations are included for 0 = 0°,5°,15 0 and 90°; 
calculations alone are shown for 0 = 30° and 55° . As additional data become available for 
0 = 30° and 55°, comparison with the predicted results of Fig. 4 should provide a 
definitive assessment of the accuracy of the present model regarding off-axis behavior. 

The point P in Fig. 3 corresponds to the set of anisotropy parameters £ 7 and £ as 
specified in Table 1. These were determined by solving (5.6)-(5.8) on a PC using the 
Levenberg-Marquardt procedure available in Mathcad 4.0. Inserting the values of 7 an d 
g from Table 1 into (5.9) gives 6 m » 85° as the fiber orientation for minimal yield stress. 

The predicted tensile curve for 9 m = 85° is virtually indistinguishable from that for 0 = 90° 
and is not shown in Fig. 4. In fact, the predicted tensile curves for the entire range 
55° < 0 < 90° show little change. Corroboration of this feature remains to be made as 
experimental data become available. 

Fig. 5 compares data and calculations for 0° tensile tests at strain rates of . 001/s 
and . 000001/s. Evidently, the strain rate dependence is captured adequately over the three 
decades of strain rate. This rests principally on the choice of the material parameter n and 

secondarily on s 0 . 

7. SUMMARY AND CONCLUSIONS 

A hydrostatic stress-dependent, transversely isotropic viscoplasticity model is 
developed by extending the well known isotropic Bodner model. The resulting model is 
tractable and retains the simplicity of the Bodner theory, particularly regarding the relative 
ease with which the material parameters are determined. It is capable of representing 
strong initial anisotropy yet is based on a single scalar state variable under the assertion 
that induced anisotropy can be ignored relative to the strong initial anisotropy. Like the 
original Bodner model, the present one is effective in representing rate-sensitive, non 
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isothermal viscoplastic responses typical of the histograms of rocket engines, eg., the 
SSME. The dependence on 7, allows for the occurrence of inelastic deformation in the 
presence of hydrostatic stress as identified by Dvorak and Rao (1976) and others. 

Natural stress states (tests) for the model are identified as longitudinal / transverse 
tension and shear and pressure- volume. Because results of shear and pressure-volume 
tests are not commonly available, an alternative characterization procedure is proposed 

using off-axis tensile tests. Conditions on £, rj and £ are examined for which a yield stress 
minimum, less than the transverse yield stress, is realized under off-axis tension. 

Although the parameters £, r\ and £ are physically based, designating the degree of 
anisotropy and hydrostatic stress dependence, empirically they play the secondary role of 
providing flexibility in correlating off-axis tensile response. A recognized limitation of the 
earlier Robinson and Miti-Kavuma (1993) formulation is that it allows matching of off- 
axis tensile data only at two orientations, e.g., 9 — 0° and 90°. Here, two additional 

orientations intermediate to 0° and 90° can be fit as well. 

Application of the model is made to a unidirectional W/Cu composite having a 
fiber volume fraction of * 9% . Characterization is based on tensile tests at 260C with 

stress / fiber axis orientation angles (Fig. 2) of 9 = O°,5 0 ,15° and 90°. Good correlation of 
data and calculations is achieved, both in terms of off-axis tensile response and strain-rate 
dependence. Predictions of tensile response at other off-axis angles, viz., 9 = 30° and 55°, 
are given. Comparison of these predictions and data as they become available will provide 
an assessment of the accuracy of the model. 

8. REFERENCES 

Allirot, D., Boehler, J.P. and Sawczuk, A., (1977), "Irreversible Deformations of 

Anisotropic Rock Under Hydrostatic Pressure", bit. J. Rock Meek Min. Sci. & 
Geomech. Abstr., 14, 77-83. 

Binienda,W.K. and Robinson, D.N., (1991), "Creep Model for Metallic Composites Based 
on Matrix Testing", J. Engng. Mech., 117 (3), 624-639 
Bodner, S.R. and Partom, Y., (1975), "Constitutive Equations for Elastic- Viscoplastic 
Strain Hardening Materials", J. Appl. Mech., 42, 385-389. 

Bodner, S.R., (1987), "Review of Unified Elastic- Viscoplastic Theory", Unified 
Equations for Creep and Plasticity, Ed. A.K. Miller, Elsevier. 

Chan, K.S., Bodner, S.R., Fossum, A.F. and Munson,D.E., (1992), "A Constitutive Model 
for Inelastic Flow and Damage Evolution in Solids under Triaxial Compression", 
Mech. o/Mat'ls, 14,1-14. 

Doraivelu, S.M.,Gegel, H.L., Malas, J.C., and Morgan,J.T., (1984), "A New Yield 
Function for Compressible P/M Materials", Int.J.Mech.Sci., 26, 527-535. 

Dvorak, G.J. and Rao, M.S.M., (1976), "Axisymmetric Plasticity Theory of Fibrous 
Composites", Int. J. Engng. Sci., 14,361-373. 

Dvorak, G.J. (1983), "Metal Matrix Composites; Plasticity and Failure", Mechanics of 
Composite Material: Recent Advances, Ed. Z. Hashin and C.T. Herakovich, 
Pergamon Press, Oxford, 73-91. 


10 



Lance, R.H. and Robinson, D.N. (1971), "A Maximum Shear Stress Theory of Plastic 
Failure of Fiber-Reinforced Materials", J. Mech. Phys. Solids, 19, 49-60. 

Robinson, D.N. and Dufly, S.F., (1990), "Continuum Deformation Theory for High 
Temperature Metallic Composites", J. Engng. Mech., 116, 832-844. 

Robinson, D.N. and Miti-Kavuma, M., (1993), "An Anisotropic Extension of Bodner's 
Model of Viscoplasticity: Application to SiC/Ti", to appear Composites Engng. 

Rogers, T.G., (1990), "Yield Criteria, Flow Rules and Hardening in Anisotropic 

Plasticity", Yielding, Damage and Failure of Anisotropic Solids, Ed. J.P. Boehler, 
Mechanical Engineering Publications, London, 53-79. 

Spencer, A.J.M., (1972), Deformation of Fibre-Reinforced Materials, Clarendon Press, 
Oxford. 


9. ACKNOWLEDGMENT 


This research was supported under NASA grant NAG-3 379. 



(C) ^ 


X 


longitudinal shear (LS) 


(a) 


(b) 


longitudinal normal (LN) 


transverse normal (TN) 



(d) 



x 


transverse shear (TS) 



hydrostatic pressure (H) 


Fig. 1 Fundamental stress states, (a) Longitudinal normal (tension) (LN), (b) Transverse 
normal (tension) (TN), (c) Longitudinal shear (LS), (d) Transverse shear (TS) and 
(e) Hydrostatic pressure (H). 



Fig. 2 Axial loading of a thin-walled tubular specimen with helical fibers. 6 denotes the 
angle between the stress and fiber directions. 
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2(^- tj) -3<f = 0 (90>) 



Fig. 3 T] vs. £, showing the limiting cases for which 6 m = 0° and 90°. The shaded area 
indicates the locus of points £, 77, £ for which is imaginary, meaning that there 
is no yield stress minimum between 0° and 90°. 



Fig. 4 Correlation of tensile data (symbols) and calculations (lines) based on parameters 
of Table 1. Data and calculations are shown for 0= 0°, 5° ,15° and 90°, calculations 
alone for $ = 30° and 55° . Strain rate is -001 / 5 and the temperature is 260C. 
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Fig. 5 Comparison of data (symbols) and calculations (lines) for 0° tensile tests at strain 
rates of 0.001/s and 0.000001/s. 


Table 1 .- Material Constants for 9% W/Cu at 260C 


n 

1.03 

m 

2.80 


L28xl0 4 (s _1 ) 

z 0 

325(MPa) 

z. 

5Sl(MPa) 

$ 

0.364 

V 

0.183 

C 

0.793 
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